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ABSTRACT 

Reasons are adduced for the continuing need for continental networks of high pre- 
cision geodetic levelling. The serious doubts which currently exist concerning the 
reliability of such networks, prompts the development of new techniques for inde- 
pendently estimating the height of mean sea level (MSL) at coastal tide gauge sites, 
Technlc^es are described for achieving these goals from satellite altimetry. 

Numerical results are obtained from the 1977 GEOS-3 altimetry data bank at 
Goddard Space Flight Center. The potential of the geoid based on the Bermuda 
calibration of the altimeter and the GM value of 3. 986 004 7 x 10^° cm^ s"^ is 
6,263,682. 8 kGalm when computed from a forty percent coverage of the oceans, 
sampled between parallels 65“S and 6 5 “N during periods representative of the 
equinoxes closest to 1976. 0, Tliis value is subject to revision by up to iO, 3kGal 
m with further orbital improvements. 

On this basis, it is estimated that the height of MSL at the Jervis Bay Datum for 
Australia is +0, 2 ^ 0. 4m. The discrepancy of zero degree between the gravity 
anomaly file for central North Amei’ica and the geoid for 1976.0 can be inter- 
preted as an estimate of the height of MSL at a Galveston Datum of +0.1 m. 

These values are in closer agreement with extrapolated oceanographic estimates 
of steric anomalies at both sites than indicated by the uncertainties given above. 
The differences in the estimates of MSLfro.m both methods are in good agreement. 

From these results, it can be concluded that all gravity data in AUSGAD 76 and 
in the Rapp gravity file for central North America refer to the geoid for 1976.0 
with uncertainties of O.lmGal. The teclinique also provides an exacting test 
of the value used for GM in the computations, 

*0n leave of absence from tlie University of New South Wales, Sydney, Australia. 
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ON THE UNIFICATION OF GEODETIC LEVELLING 
DATUMS USING SATELLITE ALTIMETRY 


1. INTRODUCTION 

Tho geodetic levelling operation Is fundamental to precise engineering worldngs 
Involving the dynamics of flow, Tlie propagation of accidental error in first 
order geodetic levelling networks is not expected to exceed ±2 -v/D^'^^mm where 
D is tho length of tho line of levelling. The systematic error accumulation Is 
expected to be about one order of magnitude smaller (Bomford 1962, p, 238), 

Tlie adjustment of such levelling networlts over continental extents should pro- 
vide data at bench marks (BMM) which are Internally consistent to better than 
+20 cm or Its equivalent. 

However, the comparison of the heights of MSL at tide gauge locations obtained 
from freely adjusted levelling networks, with those obtained from tidal analysis. 
Indicate discrepancies significantly larger than expected from the internal sta- 
tistics of the network adjustment, primarily in tho north-south direction (e, g, , 
Mather 1974b, pp, 70-71). In addition, the deduced slopes of coastal sea sur- 
face topography (SST) disagree with estimates obtained from hydrostatic 
considerations. 

In tho latter technique widely used in physical oceanography (e, g. , Lisitzin 
1974, p 72), isobaric and level surfaces are assumed to coincide at great depth 
In the oceans ( > 10 ^m). Changes dT in temperature T, dp^ in atmospheric 
pressure p_,j and dS in the salinity S of sea water (density ), produce changes 
dh in the dynamic height of the sea surface in relation to that of a standard 
column of sea water at temperature To(273‘*K), pressure p^, (1 atmosphere) and 
salinity Sq (35 parts per thousand), up to the depth of no motion h^ , according 
to the relation 



dp being the incremental change in pressure and a is the specific volume of sea 
water. 

The assumption of a level of no motion implies a coincidence between isobaric 
and level surfaces in deep oceans — a condition which is free from assumption 
only in regions exterior to the Earth's atmosphere. Hydrostatic determinations 
of dynamic sea surface heights are also subject to short period effects, such as 
frontal movements — a change of 1 mb in pressure causes a 1 cm change in sea 
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level hoight, There is also a tendency for the effect of prsesure to bo corre- 
lated with latitude (o.g, , Llsltzln 1974, p 07), It has been shown that the Irreg- 
ular lower boundary of the atmosphere causes changes In the contribution of the 
atmosphere to the hoight anomaly of ±50 cm (Anderson, et al. , 1976), The 
departures from equivalence between tsobarlc and level surfaces are expected 
to bo a maximum ns continental margins are approached. 


Global maps of dsmamlc SST prepared In this manner are not based on data col- 
loctod simultaneously. It Is recognized that the quasl-stafclonary contribution 
may contain only a fraction of the power in the spectrum of SST (e, g. , Wilson 
and Dugan 1977). 


The two factors mentioned above make some contribution to the discrepancies 
between coastal comparisons of sen surface slopes from geodetic and oceano- 
graphic methods, A third source of uncertainty Is tlie need for extrapolating 
values of the SST ( J’j) from deep oceans to coastal tide gauge sites using the 
Lagrangian equations (o.g, , Mather 1976, p 121) 


and 


X, - fX2 = 





_L ^ 

Pvf 


+ F. 


X 2 + fXi 


a?s > apa 

3X2 Pw aX2 


( 2 ) 


where (xj , X 2 ), (xj , X 2 ) and (F, , F 2 ) are components respectively of surface 
velocity, acceleration and frictional foi’ces of the ocean along the axes (X| , X 2 ) 
of a local two dimensional Cartesian coordinate system in the local horizon plane 
with the X[ axis oriented oast and the X 2 axis oriented north, f the Coriolis 
parameter 


r = 2io sin 0 , 


(3) 


0 } being the angular velocity of rotation of the Earth and 0 the latitude. 

The use of this technique assumes that current meter measurements of Xj , Xj 
are available, along with measurement of horizontal atmospheric pressure 
gradients and data for the evaluation of die frictional forces. Practical cal- 
culations are performed by assuming a non-accelerated system (i.e. , xi , X 2 = 

0). Except in abnormal conditions, Xc( < 10^ cm s"‘ and =3 0 { 0'.'2} . 

It follows diat F| , Fj must be estimated to ±10"* cm S’^ (±0. 1 mGal) if extrapola- 
tion errors are to be held below ±1 cm. Physical oceanographers have maintained 
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that extrapolations of values of over distances of up to 300 Um from deep 
oceans to coastal sites, are unlikely to Introduce errors of more than A-lOom In 
the result (o, g, , Hainon and Grelg 1972, p, 7160), 

In view of the uncertainties surrounding the osHmatos of coastal SST fi’om oceano- 
graphic considerations and die doubts cast on the validity of geodetic levelling 
(e. g. , Sturgos 1074, p, 830), It Is necessary that an Independent moans be estab- 
lished to achieve the following objectives; 

e Determination of the height of MSL at each tide gauge linked to a geo- 
detic levelling network, with a precision of at least JslOcin In the first 
Instance, 

• Definition of the universal datum level surface to which each of these 
MSL heights is referred, with an equivalent precision, 

Tliese objectives are of especial Interest to the African region where each nation 
has Its o'wn regional levelling datum, some of which have no direct access to the 
oceans and hence, the geold. Satellite altlmetrj provides, In principle, an ef- 
ficient means of achieving this objective provldeu the data Is handled with due 
regard to tlieorellcal niceties. This paper defines die basic relations which can 
provide the foundation for unifying all the world's levelling datums with a pre- 
cision equivalent to that in the radial component of nltimeter-satelllto orbit de- 
termination, On present indications, this is likely to be ±10 cm in the foresee- 
able future, In regions of adequate tracking. 


2, SATELLITE ALTIMETRY DATA IN COASTAL REGIONS 

Satellite altimetry data in coastal regions has been acquired by the radar altim- 
etor on board the GEOS-3 spaceci’aft since 1975, Tlie analysis of data in the 
Tasman and Coral Seas (Madier, et al. , 1977) in continental shelf areas off die 
east coast of Australia, provided at 2 second time intervals with the altimeter 
operating in die short pulse mode, indicates the following; 

• the sea surface appears to rise relatively steeply over the continental 
shelf slope; and 

• non-oceanic readings and hence, the transition from ocean to land, are 
clearly recognizable at the ±lm level between successive data records. 

A steep geold rlise in die region is not unexpected from the nature of die surface 
gravity field. On these figures, it can be conservatively estimated that satellite 
altimetry may provide data of quality up to 20 Ion from the coastline, especially 
when using thei ±10 cm radar altimeter planned for die SEASAT-A spacecraft, due 
for launch in mid-1978. This altimeter is expected to have a footprint of 2-12 km 
(Nagler and McCandless 1975, p, 2). 
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llio baslq data Is In the form of heights of the Instantaneous sea surface above 
the adopted roforonoo figure, llio sequence of operations to convert such data 
Into values of the heights of MSL at the regional tide gauge site is the folloivlng: 

a. Determine tlie heights of die quasi-staWonary SST In the adjacent 
continental shelf areas. This presumes that tlio goold has already boon 
defined on tlio basis of a global analysis of values of f', 

b, Extrapolate die resulting values of In the shallow continental shelf 
ocean to the coastal site using Equacion (2) and as outlined In the ensuing 
discussion. 

As fs is not greater tlian ±2m, values of should be computed from orbits 
which have a resolution of at least J:10 cm in the radial component of position. 
Values of die quasi-stationary height of MSL deduced from satellite altimetry 
should bo the average of at least one year's readings. Tlie oceanographic sur- 
veys for current velocities, atmospheric pressure gradients and frictional 
forces can only bo carried out mly on a few finite occasions, possibly just once. 
However, the continuous monitoring of local ground truth during the period of 
altimetry should provide a basis for accurate extrapolation using Equation (2), 
in most areas. 


3. BASIC RELATIONS 

The difficulties likely to be encountered in determining from satellite altim- 
eter measurements of in ocean areas have been described at length in a 
series of papers (Mather 1974a} Mather, et al, , 1976a} Mather 1978). On 
assuming the data to be of adequate precision (i, e, , ±1-5 cm in } ±3-15 uGal 
in the gravity anomaly Ag' through wavelengths of interest), the principal prob- 
lems to be overcome are the following; 

a. No complete coverage exists globally for either or Ag', Tlie pre- 
cision of oceanic gravity data is at least an order of magnitude worse 
than that of land gravity data. The former is probably subject to sys- 
tematic errors of long wavelength. It has been shown that even a 
homogeneous gravity field determination like that available for Australia 
is only adequate for SST determinations witli a precision of ±30 cm 
(Mather, et al. , 1976b). 

b. All data are measured in relation to tlie sea surface, either instantaneous 
or MSL, and not the geoid, 

c. Local MSL approximates the geoid to no better than ±2 m (Mather 1977). 
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'riireo equivalent I’olations exist between monsurod values of the gravity anomaly 
Ag', the height anomaly J*' In relation to the higher roforenoo system (Mather 
1974a, p, 91) and the son surface topography (Mather, et al. , 1976a, p, 34), 
Two of those relations are of use in determining fj through all wavelengths In 
excess of the Nyqulst limit which is defined by the nltlmotor footprint (1. o. , 
greater than about 201cm). Tliey are tlio following; 

RELATION 1 


f = (Wo - Uo J - 7(fs + koft) + Tf' - V , (4) 


where T" is die disturbing potent^^l of the solid Earth and oceans In relation to 
the higher reference system obtained by incorporating tlie gravity field model 
determined from orbital analysis (e.g. , GEM 9), with the conventional geodetic 
reference system (e, g , ,, lAG 1971) defined by the rotating oqulpotential reforenco 
ellipsoid ^vitll potential Uq ; Wq is the potential of the geoid which requires con- 
ceptual definition at the *10 cm level (Mather 1977); V is tlie potential of the 
Earth's atmosphere f;-,t the surface of measurement and 7 is for all practical 
purposes, the global mean value of normal gravity, ft Is the height of the ocean 
tide as perturbed by the Earth tide and Icq is tl\e oceanic function given at points 
^vltll latitude 0 and longitude X by 


ko(<PM = 


if (ip,X) on land 
if (<p,X) is oceanic . 


( 6 ) 


RELATION 2 

^Bcp 7N(,p/R. — 


^(V*o“Uo) “ 7 ^fsp + (\!/)[fs - fsp] /R = 

jj M,W)IN"-N"pldo. 


( 6 ) 


where Ag^ is die computable part of the pseudo-gravity anomaly Ag" on the 
Brillouin sphere (minimunri geocentric sphere containing all the Earth's topog- 
raphy), given by die relation 

Age = Ag' + 5g,, + 6Ag" + o {f Age} , (7) 
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6 g:,, being the almosplioric corrcotlon to observed gravity (Anderson, ot nl, , 
1076, p, 25), 5Ag” the ohnngo In Ag" bohvoon the Earth’s surface and Uie 
Drlllouin sphere of radius R, the subsorlpt p refer ring to values at the point of 
computation, do in Equation (0) rotors to Uio element of surface area at an 
angular distance fp from tlio point of computation, 

'Pile gravity anomaly Ag' on the higher reference system can bo related to llio 
conventional gravity anomaly Ag by Uie relation 

Aij' = Ag - 67, (8) 


where 67 is defined by the coefficients C'^mn Earth gravity field model 
like GEM 9 (0. g, , Loroh, et al, , 1977) by a relation of the form (Mather 1974a, 
p. 96) 

n' a 2 

g7 = 7 ^ (n-1) X] ^ ^aiiin + e , (9) 

n=0 m=0 a=I 

S^nin being surface spherical harmonic functions of degree n and order m, given 
by 


Slum “ P„i„(sin<^) cos tv.N: S 2 ,„n = P„„,(siii 0) sin mX ; (10) 

n' being approximately 20, 

Tlie quantity Nj.' is related to the height anomaly -■ l.e, , tlio height of the sea 
surface above the reference surface In ocean areas - by the relation (Mather, 
et al, , 1976a, p. 29) 


-1(V-6T"), 


(11) 


6 T" being the change in T" between the Earth's surface ard the Brillouin sphere. 
Tlie surface Integrals in Equation (6) apply on the lirilloulu sphere which is the 
smallest sphere on which the ortliogonal properties of surface spherical har- 
monics apply without approximation, to the potential T" and the pseudo-gravity 
anomaly Ag". 

The kernel Mj(0) of the surface Integrals in Equation (6) is given by 

00 I 

Mj(0) = ^ n(2n + 1) Pno(cos 0) = - -j cosec3*/20 - 3 cos0 if 0 0 . (12) 

ll'=2 ^ 
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Equations (4) and (0) take Into account tho fact that data is recorded In relation 
to the sea surface and not the geold. Equation (5) can bo applied to tho instan- 
taneous value of but Equation (0) has to bo evaluated from global stationary 
fields of i*' and Ag'. It is Ujorofore assumed Uiat tho effect of tides has boon 
removed from Hie data pr*'/r to use in numerical evaluations. For methods In 
handling the lidos in altimetry data, see {Zotlor and Maul 1971{ Brotrogor 197G; 
Mather 1077, See, G). In the case of data on land, tlio value of at tho element 
of sur face area da refers to tlio height of MSL above tho goold at Uio regional 
levelling datum. 

On considering tlio shortoomings listed at tho commencement of this sootion, Uio 
most favorable procedure for dotonnlnlng fj appears to bo tlio following (Mather, 
ot nl. , lOTGa*, Mather 1978), 'ITio spectrum of non-tlrlal quasi-stationary SST 
( fs ) can bo considered to bo constituted as follows 

fs = fss fse . (13) 


where are components with wavelengths longer than tliat (fi) In tlio Earth's 
gravity field which perturb altlmoter-aatollite orbits above tlio noise level of 
tlio tracking. It Is estimated that fi = o (lO^ Ion} for 800km altitudes whore 
the satellite is tracked with ±10 cm all-weather systems from a global networli., 
refers to all conti’lbutlons with shorter wavelength. 

The contributions fjg to the SST can be determined from the following 
equations} 


EQUATION 1 


= (i)"t 

' 11“ 0 \ / ni«o a*l 


7r+(Wo-Uo)-7(fs+W- V. 


(14) 


where M is tlie mass of tho Earth with atmosphere, G tho gravitational constant 
and are harmonic coefficients of degree n and order m, Tho first equality 

applies In the space at and exterior to tho surface of measurement while the 
second applies at tho surface itself. All coordinates (R, ^ , \) are geocentric 
spherical coordinates. 
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rQUATtON 2 



where 

Agd = Ag' + figg - i^gfd + 2 G 0 + o {f^Ag'} , (16) 

being the lefleetlon of the verfiGal at the point where gravity (ig) is measured 
and is given by 

G 0 = m ^ f ^ 3f sin? 0 , (17) 

where m is defined in Equation (46). flie first equality at 15 provides &e defin- 
ition of Ag" , the second applies in the space exterior to &e surface of measure^ 
ment while the tMrd applies at the latter. 

Equations (14) and (15) can be used to determine fj with wavelengths greater thsm £ 
When the coeffieients Ganji, iii Equation (9) are known to the equivalent of 
o {iO.Oi^O, OS ldjalm} . In practice, the satellite deterniined gravid field 
model is already incorporated in the ht^er reference System and hence reflected 
in the values cOmputeu for Ag' and f Thus 

Canm = ® o {ir»} for n<n . (18) 

The resulting observation equations take the form 


Equation 14'- for values of f' from satellite altimetry 


Vf' ^ 


V 

-y 


kgft ■* 

■i(Wo»Ui 





\ 


Order of ( m) ±4 

Magnitude 

44 

known 

5 tty? 
known 

i2 

local 

value 

±1 

local 

value 

*1 

zero 

degree 

Range of (co) ^ ,, 

Wavelength 

0 < OJ <o® 

0 <00 <“ 

oo>£ 


= oq 
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and 


Equation 15 - for values of Ag' in land and continental 

shelf areas 



VAg' 

= “ 

^ (1 “ " 

27 2 

(Wo-1 



4- 


•I i 

4- 

Order of (mGal) 
Magnitude 

±10 

.MO 

known 

±Yi 

unknown 

constant 

±Vl ±V4 

local local 

value value 

zero 

degree 

Range of (w) 
Wavelength 

w<2 

0<CO<oo 

O)~oo 

co>£ 

= oo 


Tho baeiK for the recovery of the SST by the use of these equations under con- 
ditions of unfavorable signal- to -noise is the band-limited nature of the signal 
being recovered, as can be seen from the wavelength ranges listed above. 
Equation (20) also constitutes a basis for recovery of the MSL height at the 
regional elevation datum ( ) as discussed in Section 4, 

4. CONDITIONS INFLUENCmG DETERMINATIONS OF THE HEIGHTS OF 
MSL AT COASTAL SITES 

The numerical value of depends on the basis adopted for the definition of 
the geoid. Two possibilities exist (Mather 1977): 

a. Adopt an oceanic definition for the geoid 

The geoid in such a case, is the level surface corresponding to MSL in 
ocean areas, such a definition can be realized by representing the SST 
by the relation 

n' n 2 

fs “ S 13 ^sanm ®anm fss » 

n=I 'm=0 «»1 

on using Equations (5) and (13), SQ^n, being given by Equation (10) while 
fsanm surface spherical harmonic coeffiolents of degree n and order 
m, n' having the same significance as in Equation (9), For a more 
comprehensive discussion of this problem, see (Mather, etal, , 1978, 
Sec. 6). On defining by the relation 

+ (5T" - V)/7 (22) 


original page la 

OP POOR QUALIT*) 
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observation equations of the form 

Vf' = - ko(0,X) Samn " ‘ l<o(<iA)) f«d + “ " ^o) (23) 

n-i \ / y 

nre sot up, with the apbericnl harmorio eorios only reprosontlng In 
oconnlo nroae. In land areas, kg = 0 and the effect of SST is represented 
by the value i'sjj at the regional levelling datum, 

Ulio observation equation for grravlly anomalies takes the form 

''Ais* “■ kgf0,\) i^sanin^anm “ (Wg - Ug) (24) 

on using Equation (21) In Equation (20). The value of Wo - Uo can bo 
considered to be a Imown quantity. Its magnitude Is dependent on the 
definition adopted for tlie geold, as discussed In Sec. 6,1. 

Under such conditions, and as the spherical harmonic modelling and the 
discrete values are mutually exclusive In Equations (23) and (24), It eon 
bo argued that fsj can bo determined solely from the amlysls of data re- 
latcd to the regional levelling datum, 

b. Adopt a boundary value problem definition for the geold 

In tills case, fs is defined by tlio equation 

n' n 2 

fsanm ^anm ■*“ i*s» ( 26 ) 

n“l n)*0 ot=I 

instead of Equation (21), 'ITie value of f generated on land from such a 
model should bo constant for locations on the same regional levelling 
datum, equal to tlie ^sd value for the datum. It would not be unreason- 
able to expect aliasing effects in coastal areas (ibid. ) and the portents 
aro not favorable for obtaining estimates of fjd datums with extents 
smaller tlian B^Ion^, 

. -7 , • . . 

It con be concluded' that the adoption of a specific definition for the geold fixes 
the value of tlie zero degree term (Wo - Uq). This can only be obtained from the 
global analysis of values of (and Ag', If the second definition Is adopted for 
tlie geoid). It also follows tliat the analysis of gravity anomalies on regional 
geodetic datums can provide a basis for the determination of the value of fjd 
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at Uio regional lovolllng datum If the area covered by the datum ts larger hi ex- 
tent km^) than the highest full harmonic In the higher reference model which is 
free from orror. The only role played by the sateUlto altimetry is in defining 
(Wo - Uo). 

In using Equation (24) for land areas (ko “0), the density of gfravlty data points 
in well reprosontod areas is 1 per 10^ km while the number of high vireoision 
position fixes at which Equation (23) can bo used is more likely to bo 1 every 
10® lon^ . The extent of high frequency noise in the more probable gravity 
anomaly observation equations can bo reduced by fonnlng the latter using larger 
area moons, say two degree area moons. As the values of ouch area moans are 
strongly correlated with position (e, g. , Mather 1975, p, 77; Mather, et al, , 
1976b, p. 78) and in view of the adverse signal-to-nolse illustrated In Equation 
(20), it is prudent to model those variations in V/^g' , which are two orders of 
magnitude larger than tlie contribution of the term containing fsd • Any two 
dimensional model should siifflce for the task, assuming that the gravity data is 
evenly distributed about the datum. Thus Equation (24) can bo written for land 
areas on the same levelling datum in the form 

Vfis' = Asa - + E £ £ Ci'„„ + 2 (W„ - U„) , (26) 

where are harmonic coefficients of the Fourier fimotlona Fanm * defined 
by 

Finn, “ cos(n + m AX); F2nm “ sln(n A0 + m AX), (27) 


both n and m not being equal to zero simultaneously. A^, AX In Equation (27) 
are differences of geocentric surface coordinates from some convenient point of 
referencvs in the region. The moat important wavelengths which need to bo 
modelled in order that the resulting value of is not aliased, are die 
following; 

a. Those equal to 4 times the smallest dimension (d) of the region served 
by the datum arising from errors in the assumption described by 
Equation (18); 

b. Those equal to twice d, due to residual errors in the gravity and 
levelling networks; and 

c. The error of assumption at 18. 
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Tho quality of tho determination will depend on 

o the extent of the area served by tho levelling datum, represented In the 
solution; and 

* wliotlier nil wavelengths shorter than H have been sampled In the 
determination. 

Anolher factor Influenoing the determination Is tho precision with whloh tho 
regional gravity datum has boon ostabllshod, Tho required precision for a 
±10 cm determination is ±30 pGal - a precision which can be achieved with 
modem transportable absolute gravlmetors. 

It also follows that the quality of the determination will diminish as a function of 
tlio shortfall below fi^km^ of the area served by the regional levelling datum, 

'riie use of Equation (24) cannot be expected to give stable results if the area 
sampled is loss tlian ‘AC^km * , oven If tho modification at Equation (26) were 
used. In such cases, It becomes necessary to resort to Equation (2) for ex- 
trapolating fs from oceanic determinations using satellite altimetry as dis- 
cussed in Section 2. 

A determination of on a global basis using Equations (19) and (20) will only 
define the contributions fjj in Equation (13), There Is no reason to believe 
that tlie magnitude of fjs is any Smaller Uian that of , A preferred pro- 
cedure for establishing at tide gauges serving as levelling datums for regions 
smaller than kiC^lcm^ or at those not connected to levelling networks Is the 
following: 

iVflgr 7 

Define fsj at a set of locations about 20 Ion offshore using Equations (19) 
and (20), 

S/age 2 

Having accomplished Stage 1, define fjs using Equation (6), 80% of the 
contribution made by the surface integral Is expected from regions In the 
range 0?1< ^ < 5® If tho higher reference model Is used (Coleman and 
Mather 1976). The balance 20% comes largely from tho Innermost zone. 

As (Nc - Ngp ) In Equation (6) can take botli positive and negative values in 
tills region, the possibility exists tliat the band-limiting constraints placed 
by the finite footprint of the altimeter, results in the loss of high frequency 
signal. The significance of this depends on the amplitude of very short 
wavelength contributions ( < 20 km) to I's. This effect can be disregarded 
If less than ±5 cm, noting that the quantity required Is the annual average SST, 
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Stage 3 


Use current meter data, atmospheric data and modola of frictional forces to ex- 
trapolate obtained at oceanic sites within 20 km of the coastal site in Stages 
1 and 2, to the latter using Equation (2). Tho resulting value should agree ^vlth 
that obtained from processing all gravity anomaly and height anomaly values re- 
lated to tho datum, provided adequate coverage existed for tho purpose. 

Solutions for in Stage 2 are also influenced by the occurrence of tho term 
within a surface integral, as discussed in (Mather 1977, Sec. 7). It la therefore 
necessary to solve for fgj as a regional field using models of the type described 
in Equations (2G) and (27). This, in turn, requires that Ag' be defined as a 
regional field in this basic oceanic area wiin a precision of ±0. OSmGal thi’ough 
wavelengths of interest. In theory, the minimum wavelength in fgj which can 
be recovered, is governed by tlie Nyquist limit which Is a function of the satellite 
altimeter footprint. In practice, it is desirable to reduce the frequency range 
in fss in view of the difficulty of defining Ag' to the required precision, free 
from tho effects of variations In f' witlx periods shorter than that implied in tho 
quasi-statlonary concept (i, e. , less than a year). It would be most helpful if it 
were established that quasl-statlonary SST over fills period had a power spectrum 
to which tho contributions of wavelengths below some lower limit (say, 10^ km) 
were less than 20 cm^ , 

Another problem In the evaluation of at Stage 2 Is the necessity to define Nc 
(Equation (11)) as a continuous field within 600 Itm of the point of computation. 

This would call for the determination of to ±10 cm in all land areas which 
fall within this region. The only means of achieving this objective at tho present 
time, is by collocating transportable laser traeWng systems at points in the first 
order geodetic levelling network in the area. Ideally, (and hence Nj.' ) should 
be defined at points on a 10^ km grid. The cost of such an operation la 
prohibitive. 


5. PRACTICAL CONSIDERATIONS 

It has been shown in Sections 3 and 4 that satellite altimetry has the potential 
tos 

a. select a particular level surface of the Earth's gravity field as the geoid 
at the ±1-5 cm level by defining a magnitude for (Wo - Uq ); and 

b. define the height of MSL above the geoid so selected at coastal sites in 
conjunction with gravity, oceanographic and geodetic surveys, 
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lie pFecision requirements to be met by the varloue (ypes of data are tte 
following; 

a. Sateliite altimetry data . 

Sea surface helots' must be bomput^ from o Altai epbemerides 
with a FesolUtlon of cm in the radial cOmp^ent of position* 

b. Global gravity field model to degree n' 

• * • I ' 

I%e global gravl^ field model consistent with the above oAltal ephem^ 
erldes to the noise level of the traeMng should be suoh ^t When used 
to define the higher reference tuodol, ^ resulting values of T" have 
no terms with wavelengQis greater thkn f (equivalent to n' ), £ being 
estimated to be o { 10^ ♦ > 

c. Regional levelling surveys wifhiti iOQ'km 0 / the computatibnal area 

The precision of such control [mf Veys used In ^ control of gravity 
anomaly and height anomaly cbmputatlpns, should cover coottnentai 
extents. 

d. Regional gravity control networks 

The precision required on land is ±30 nGal In ail control networks with 
wavelengths greater than £. 

e. Oceanographic surveys for extrapolation to coastal sites 

The precision required for cur rent velocity measurements Is 1 cms~i 
and for frictional forCes 1>0~^ cm S~^ If a ±1 Cin reSolUtlott is to be ob^ 
talned in the extrapolation of fg . The values required are tte average 
for the sanipllng period. 

Ihe desired precisions are not ctmrently available for data at (a)* (b), mid (d). 
in the case of data at (e) . the favorable indioatimis from die Internal statistics 
of level net adjuspnents need to be temi^red by the doubts Implied In estimates 
of coasM SST eSpeci^ly in the merldion^ dlrecdon^ as obtained from levelling/ 
tide gauge comparisons. ^ 

An attempt to improve the gra^diy data bank for Australia for use In SST deter ^ 
mlnations resitlted in a gravity ahomaly representatloh which was assessed as 
being sufficient only for determinations to ±30 cm In . (Ma^r» et al. , 

i 976b) . None of the other data btmks a A likely to be of better quality. 



The quality of global gravity field models currently available, has boon tlio sub- 
ject of close scrutiny, especially vdth the advent of satellite altimetry data of 
quality from the GEOS-3 mission. Recent studies of GEOS-3 altimeter data off 
eastern Australia show discrepancies between geoldal and sea surface models 
which appear to have amplitudes of up to 5 m o;nd wavelengths of 2000 km (Mather, 
ot al, , 1977, p. 36), On the other hand, it is possible to manipulate GEOS-3 
altimetry data so that the dominant features of the global quasi-stationary sea 
surface topograpliy are recovered witli a precision estimated at ±0 cm (Mather, 
et al, , 1978, Sec, 9). It was tht refore decided to use the techniques described 
above to: 


• define a geold consistent witli the 1977 GEOS-3 altimeter data bank; and 

• establish, if possible, the height of the datum level surface Implicit In 
the current gravity data banks for Australia and the United States, 

Wlille some doubt exists about the practical significance of the latter results, 
the computations would highlight the nature of the numerical problems encoun- 
tered in the evaluation. Those computations are described in Section 6 and the 
results discussed In Section 7, 


6. NUMERICAL RESULTS PROM GEOS-3 DATA 
6.1 The Geold for Epoch 1976,0 


See Section 4, As abnormal conditions may prevail in coastal areas, it Is pre- 
ferable to select the datum level surface on the basis of data sampled in ocean 
areas alone. 

The Data Set 

Tlie GEOS-3 altimetry used in the definition of the geoid i&r epoch 1976. 0 were 
the total data available in the 1977 data bank at Goddard for the periods 1 Sep- 
tember to 31 October 1975 and 1 March to 30 April 1976 (the Equinox Data Set 
for 1976, 0). 350 passes of data were recorded during the first period and 284 
during the second. The total distribution of data is shown in Figure 1. This 
data was used to obtain a geometrical model of the sea surface consistent with 
the best available orbits as described in (Mather, et al. , 1978, Sec, 5), The 
resulting model is based on data which is minimally affected by the seasons. 

The representation obtained, however, is less than desirable due to the irregular 
data acquisition. 


15 


Tho Computational Procedure 

The instantaneous position X| of tho son surface at epoch (t = t) is obtained from 
tho Earth space satollito coordinates Xjj (t) on a tliroe dimensional Cartesian 
coordinate system X{ related to the geoconter, tlio CIO polo and the meridian of 
zero longitude (Greenwich), and the gravitationally stabilized altimeter range 
lii(t) to tile sea surface, by the relation 

Xj(t) - Xis(t) - h(t)Ki(t) + o {l0-8ii} , (28) 

where je{(t) are tho direction cosines of the normal which passes through the 
satellite position. The coordinates Xf are easily converted to geocentric 
spherical coordinates (Rq, 0, X) by the well Icnown relations 

1^0 = ''"*5 ^ = tim-»(X 2 /X,); 0 = tan- 1 [XjAXj +^f]. (29) 


If the coefficients of the Earth’s gravity field model are "free from error" to 
degree n' , the geopotential at the sea surface Wss can be modelled using the 
harmonic coefficients (n < n' ) using the relation 


Gm 


- R, 


u / \ n n 

E l ^ l" 
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+ }4RoCos^0£o2 , 


(30) 


111 “0 a= I 


where are defined by Equation (10), and a Is the Earth's equatorial radius, 
Tlie m,cdel for tho potential Vj of the atmosphere extei’ior to it is given by 



^sanm ®anm > 


(31) 


where Rg is the radius of the minimum geocentric sphere enclosing the Earth's 
atmosphere. Downward continuation through the atmosphere is only possible 
when considering the potential We of the solid Earth and oceans, given b3' 


= (ly z E c.s, 

^ n=D \' / m=0 a= l 
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an hi > 
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where 
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Tho goopotontial Wjs fit tho soa surface is given by 

Wss = CWe)„ + Vg, (34) 

where Vg is the potential of tho atmosphere as evaluated at tlie sea surface. As 
tlio atmospheric potential Is dominated by the term of zero degree (Anderson, 
et al, , 1976, p. 33) and also contains harmonics of degree 1, considerable com- 
puter economy is achieved by representing Vg by a surface harmonic model 

OQ tl 2 

Vg Vgattjn ^anm o {0.02kGal m} . (36) 

n°o m“0 


Thus the low degree harmonic estimate (n < n' ) Wss of the geopotential W at tho 
sea surface can bo written as 


w V 

° n=0 
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E 

n- 1 


c' 

an m 


^anm 




(36) 


as evaluated at (Rq > ^ > X) defined by Equation (29), where 

r = r' + sr 

'^anm '-'an in ^ ^'“anm* 


(37) 


The coefficients Ck,,,,, in Equation (37) are satellite-determined harmonic co- 
efficients of die type embodied in the model GEM 9 (Lerch, et al, , 1977), while 
5 are the corrections required in downward continuing tho satellite de- 

termined geopotential for deteimining Wss ut the sea surface , being defined by 
the relations 


Ro 


Table 1 sets out die corrections computed from Anderson's evaluation of 

the atmospheric potential atbotli satellite altitudes and the Earth's surface 
(Anderson 1976), Tho effect of the corrections to Wss per coefficient never 
exceed UGalmm for decrees up to (6,5). This is due to the gravitational effect 
of the long wave componi-nts ( > 2000 km) contribute 98 percent of the strength 
of signal which vary less tlian ±5kGalcm over the surface of the oceans (ibid. , 
pp, 209-210). Consequently, the differential effect is insignificant through low 
degree terms. 


17 


Under Uroso olrouniBtancos, W^s can be oomputed from tb© coofflclonts Cctiiin 
Instead of In Equation (30). 

Evaluation of Wq 

Tlie 034 passes of GEOS-3 altimetry between latitudes -GB'’N and C5®N recorded 
In Soptembor-October 1976 and March- April 1970 provided a 30, 8 percent rep- 
resontation of the 33,902 oqul-angular 1“ x 1® squares olassiflod as bceanlc in 
this study. 'Hie rms residual representing variations vrtthln a square was 
M. 4m, Tlio resulting value of Wq obtained was 

Wo = 6,263,682.76kGahn, (39) 

based on the Bermuda calibration of the GEOS-3 altimeter, and 

GM = 3.986 004 7 x 102°cm3s■^ (40) 

consistent witli tlio velocity of light c taldng the value 

c = 2.997 924 58 x lO^^cins"'. (41) 


This value has an estimated uncertainty of ±0.4m. For details on how the data 
was processed, see (Mather, et al, , 1978, Sec, 3). (A summary of results is 
given in Table 4. ) 

0.2 Tlie Reference System Used 

Tlie system of reference is defined by the gravity field model GEM 9 whoso 
coefficient 


C 20 = -1.082 627 (6) X lO’^ 


(42) 


is consistent with the value of c at (41). On adopting 

a = 6,378,140.00 m, (43) 


it follows til at the potential Uq on tlie surface of die rotating equipotential ellip- 
soid of reference is (e. g. , Mather 1971, p. 83) 


U„ 


GM 

a 


a 
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(44) 
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( 46 ) 


being related to C 20 by the relation (o,g, , 

silica r 



Mather 1978, App,) 
2 111 siiia 
15 q2(a)^ ’ 


where a = oos“* 1 1 - f ) , 


in = u^co^/GM, 


and q 2 (a) is given by 


qz(«) 



[a (3 cot^oi + 1) -* 3 cot ot) . 


m 


(47) 


The referonoe quoted sots out simplified procedures for solving Equation (46), 
Uo for the system of reference adopted In tlie present series of calculations Is 

Uo = 6,263,682.67(6)kGalm . (48) 


6, 3 The Computation of the Gravity Anomaly 

Tlie correct procedure for preparing gravity data for high precision compu- 
tations is described In (Mather, et al. , 1976b). The data maintained in most 
gravity data banks (e, g, , Rapp 1977) is in the form of free air anomalies Agf , 
on some system of reference, usually Geodetic Reference System 1967 (GRS 67), 
computed from the formula 

Agr = g - To + 0.3086 h(”0 (49) 


where Tq is normal gravity computed for tlie oqulpotentlal ellipsoidal model 
using the formula 

To = Te(* +/3siii2^g +/32sin2 20g) + 0 jo.OSniGalj , ( 50 ) 


where 


Tc = 978,031.675 niGuI 
(i = 5.30254 X 10-3j jSj - -5.862 x 10'^^ . 


(51) 
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As pointed out in (Ibtd,). resolution to ±4,uGnl onn bo obtained by using Uio 
formula 


y^i ° y^O +|3'siii20g +/32sin‘^g), 


( 62 ) 


where 


Tc ■= 978,031.678 mGal 
|3’ = 5.27893 X lO*’; |3^ = 2.346 x 10*'^ , 


(63) 


Tlio gravity anomaly Ag has to bo re-computed using the formula 


Ag - Agf - 0.3086hC»') 


2 AW 
u 



> 


(64) 


whore is defined by Equation (17) and AW is related to Increments of geodetic 
levelling dz by the relation 

g dr . (55) 

1)1 



D a t ti 


For a description of those calculations for the Australian gravity data bank, see 
(ibid, , p, 68), The gravity data for the United States was In the fonn of free air 
anomalies originally computed by the Defense Mapping Agency Aerospace Center 
(DMAAC), 'Hie normal gravity was computed using closed formulae (e.g. , 
Mather 1971, p. 88). Conversion to gravity anomalies was obtained by using 
Equation (55) on the gravity and elevation data banks for the Central North 
America and using the resulting set of geopotential differences related to the 
Galveston tide gauge. 

All gravity anomalies wore finally referred to the higher reference model defined 
by GEM 9 and the constants defined by Equations (41) and (43). 


7. COMPUTATION OP DATOM LEVEL SURFACE DISPLACEMENTS FROM 
THE GEOS-3 GEOID FOR THE EPOCH 1976.0 

7.1 The Jervis Bay Datum Level Surface 

All Australian gravity data is related to a freely adjusted Australian Levelling 
Survey of 1970 and referred to the Jervis Bay Datum at (<;^ = 35. l^S, X 150. 7“E) 
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as dosorlbocl In (Mathei’, otni. , 1970b, Soo, 2), In a first stage, the gravity 
anomalies Ag wore converted to th pseudo-anomnllos Agj using Equation (10), 
In preferred clroumstancos, the liigher roferonco model used should bo free 
from error. Assuming that is required to ±10 cm, it follows Umt Equation 
(18) must also hold for all coofficlonls Included in the higher roferonoo model, 
assuming die number of coefficients to bo about 400, A solution procedure based 
on Equation (26) will bo subject to considerable aliasing of the value of ?sd If die 
errors in the higher reference model with wavolongUis greatej tlian tlie shortest 
dimension (d) of the area served by the datum, were not modelled in die compu- 
tations. It is ostimated that the error In the GEM 8 coofflolonts to (4, 4) on 
models at the surface of the Earth is ±1.4 x 10‘® (Leroh, ot al, , 1977, p, 52), 
equivalent to approximately ±9 UGal om in t'' , Tlieso ostimated errors increase 
rapidly widi increase of n to around ±60kGalcm for degree 20, 

The value of can in principle, bo obtained by the analysis of either the 1* x 
1*, 2® X 2“ or 5® X 6® data banks. Tlie results obtoinod are influenced by die 
following factors! 

• Tlie signal-to-nolse . fsd is not larger than ±2 m while the variability 
of the data Increases V'ith decrease of square size (Taiilo 2, Row 3). 

• Departures from the assumption that the gravity field model is orro ?' 
free . Tlie existence of a largo non-zero value for the regional moan 

( Ag^) of Agj over Australia emphasizes the need for Fourier model- 
ling the long wavelongdis errors in the gravity field. The large positive 
values of Ag j for Australia (Table 2, Row 2) indicate die net high of 
surface gravity in die region, Tliose values are highly oorrolated widi 
position slioivlng net highs in the oast and west of the continent ivith a 
band of lows In the center (e.g. , Madier, et al. , 1976b, p. 78; Lorch, 
ot al, , 1977, p. 71). This type of effect has a wavelength tivo-thlrds 
that of (ho east-west dimension of tiio continent and should bo modelled 
when using Equation (26), 

• Errors in die area means . Those arise primarily due to inadequate 
sampling. 

It was therefore decided to model the following wavelengths in the Fourier series 
when effecting a solution: 
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Tlio values adopted for d In the Australian calculations wore d^ « 30* in latitude 
and d;^ “ 45* in longitude. 

Tn viow oC Uio unfavorable signal-to-nolso, it was necessary to constrain the 
solution to on a pilori assessment of the mag;nltudos of the oorroctlons, Vor 
example, the toim tSsu /R in Equation (2G) will not exceed ±0, 3mGal while the 
coofficionts C'a,,,,! should on the average, not bo significantly larger than Ag,j/N, 
where N is tlio total number of harmonics modelled, ConsequouMy, the solutions 
sho\vn in Tables 2 to 4 were obtained by minimizing 

n' n 

4' = 13 I] Wc|(Ca,„n)2 . (6G) 

l«t i**! 

whore 

I 100 if a« 1, n « in ° 0 

W| " l/cos0i: Wci = (57) 


The solutlono obtained for Australia using AUSGAD 76, GEM 9 and the GEOS-3 
alllniuEer-detormincd geold for 1970, 0 are sot out in Tables 2 and 6, Tlio pre- 
ferred result is obtained using fully repx’osonted 5* x 6“ area moans as the area 
moans ar'ii probably more reliable, being loss affected by Irregularities in 
gravity field sampling. Tlie number of observation equations is limited, reducing 
to 15 if only fully represented squares (1. o, , only 6“ x 5* squares based on 25 
1* X 1® values) were considered (Table 5), 

On this basis, the preferred value for the height of MSL at the Jervis Bay Datum 
is 


(fsd)jcrvis Bay “ - 


(58) 


The equivalent value as extrapolated from the deep oceans using oceanographic 
data is +0.1 ± 0,2m (Mather, et al. , 1978, Figure 1), noting that a zero degree 
effect of +1. 14m has been eliminated. TJie figure at (58) is referred to epoch 
1968. 0. Tlie variation of the height of MSL with time at Sydney Is estimated at 
+1 mm per year. Thus there is less than 1 om discrepancy introduced Into the 
result due to the non-coincidence of epochs of the levelling and the altimetry. 
'Hie error In the datum for the Australian gravity is estimated at ±0. 06 mGal 
(Mather, etal. , 1976b, p. 79), introducing an uncertainty of 0.18 m in the re- 
sult at (58). For estimates of other sources of error, see comments on the 
result at (59), 
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7,2 Eglimating U 10 Effeots of Zero Degree In the Gravity Data Bank for Central 
North Amoiica 


Tlie region covered by this study was Uie North American continent bounded by 
the parallels 28“N and 50®N. Tliis Included a small part of Mexico and the souUi- 
eastern part of Canada, Gravity values on the North American continent are, as 
best as possible, referred to the International Gravity Standardization Network 
(IGSN 71) (Morelli, et al, , 1971), The basic network was assembled by the 
Defense Mapping Agency Aerospace Center. Tt would be difficult to assess, 
without a major ro’-examinalion of the data, whetlier the pattern of errors in the 
United States Levelling Network are reflected in the resulting 1® x 1“ free air 
anomaly data bank compiled by Bapp (1977). The data used in this study had been 
rounded off to the nearest mGal, Its characteristics are summarized in Table 
3. Parts of the Canadian gravity data bank were also included in this study. 

Tile same comments made about die elevations of gravity stations in the United 
-States apply to those in Canada, there being a variable systematic difference 
between common jimction points of the two levelling systems which is about +10 
cm on the average (Lachapelle 1978), This has not been considered significant 
in the present study which is of an exploratory nature. 

It is therefore not clear that the analysis of the gravity anomaly data bank for 
central North America, prepared as described in Section 6. 3, ivill contain any 
information on the height of MSL at the datum level surface for the region, as 
implied in the computation of free air anomalies. As a starting point, it was 
decided to adopt the suggestion that tide gauge at Galveston be adopted as a suit- 
able datum for this levelling network (Holdahl 1978), Geopotential differences 
were computed using 1® x 1® mean square elevation and gravity data banks in 
relation to the value in the 1® x 1® square (0 = 29.5®N; X = 2G1,6®E) containing 
this site, 

Tliese data banl?s and the resulting geopotential network were used in Equation 
(54) to produce a gravity anomaly data bank for central North America as de- 
fined above. The characteristics of the data used in the analysis are shown in 
Table 3. As mentioned in the previous seclion, the discrepancy between the 
GEM 9 model and the surface gravity data, as embodied in the value of Agj for 
the region is five times smaller than that for Australia (Tables 2 and 3, Row 2). 
This is probably a reflection of the better tracking coverage available in the 
North American area when compiling the GEM 9 model. 

If it were assumed that all the gravity data in the Nortli American study were 

« based on a regional standardization network of the same quality as 
IGSN 71; and 
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• convortod to gi'avity anomallos based on a network of elevations sub- 
stantially controlled by the freely adjusted regional levelling network, 


it can be said that the height of MSL at Galveston is given by 


(fstpGiilvcslon 


“ +0rM m . 


(B9) 


It Is only possible to obtain a very approximate oooanographical value for fsd at 
Galveston as +0,1 i- 0,3 m (Levitus and Dort 1977, p. 1283), allowing for the 
zero degree effect. The sources of uncertainty in the result at (59), provided 
the above assumptions were valid, are the following: 

• ±20 cm due to error's in the gravity standardization network. This 

figure is a guess, compatible with the more carefully assessed figure 
for the Australian national network, quoted in Section 7, 1, 

e ±12 cm due to aliasing as a result of using too few coefficients In the 
Fourier modelling - an inevitable consequence when using area means 
for improving the si gnal-to -noise, 

s Tlio value of Wo obtained in Section G, 1 was not based on a full coverage 
of die oceans between 66®S and 65®N, As shown in Table 4, die result 
may require revision by up to ±30 cm as further orbital refinements are 
made. 

It is not unreasonable to conclude that the values of fsd given in Equations (58) 
and (59) have uncertainties at the ±0. 4m level, Tim level of agreement obtained 
^vith oceanographic values is much better, being about one-fourth this value (l,e, , 
±0, 1 m). 


8. DISCUSSION OF RESULTS 

Tlie results presented in Section 7, 1 and 7, 2 are based on the following data; 

• A geoid for epoch 1976. 0 based on data in die 1977 GEOS-3 altimeter 
data bank. This data base is being added to and in die process of further 
revision. It is not expected that the value of Wq given in Section 6. 1 
ivill change by more than ±0, 3 kGalm when the representation increases 
from the 39,6 percent coverage used in the px-esent study and when 
refinement is complete. 
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• Tlio gravity anomaly data bank for Australia spoclally prepared for sea 
surface topography determinations (AUSGAD 76). 

• Tlio 1® X 1“ free air anomaly data set for central North America origin- 
ally compiled by the DMA /'verospace Center and provided by Rapp in 
the form of values rounded off to the nearest mGal, 

• The GEM 9 gravity field model. 

Tlie last named data sot is not critically Involved in the determination though 
wealmesses In die model cause additional signal-to-noise problems, 

Tlie following observations can be made on the results presented in Sections 6 
and 7 s 

a. The analysis of tlie data for Australia (Tables 2 and 6) indicate die ex- 
tent of the aliasing influence of 5® x 5® area moans which were not based 
on a full representation of surface gravity data (i, e, , twenty five 1® x 
1® values). Restriction of the analysis to fully represented areas re- 
duces the ratio of unlmoivns to observation equations. This is offset by 
the reduction of noise in the observational data and results in an im- 
proved solution. Stability of solution is enhanced by restricting the 
Foiuier modelling to tlie same range of longitude per parallel sampled. 

b. The results for Australia indicate diat the use of this teclmique in 
regions not providing heavy tracking coverage for the development of 
the satellite determined gravity field model, will produce conditions 
where fsd to bo determined in the presence of adverse levels of 
noise. Subsequent computational instability can be avoided by studying 
the nature of the distribution of Ag^ over the region before the selection 
of wavelengths for Fourier modelling, 

0 . Tlie results given in this paper for the MSL datum at Galveston are based 
on the assumption that the gravity anomaly data bank for central North 
America was based on the geodetic levelling. There is no assurance 
that this is the case. It is most desirable that this experiment be re- 
peated with a gravity data set whoso elevations are Im'^wn to bo related 
to die continental levelling network based on the Galveston Datum Level 
Surface. 

d. The results presented in this analysis establish the potential of this 
method for defining the height of MSL at die regional levelling datum 
serving areas larger dian the square of the minimum wavelength in the 
satellite determined gravity field model. Ideally, the model should 
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be free from error tlirough these wavelengths. However, slightly de- 
graded results can bo obtained even If this condition Is not satisfied, ns 
seen from the results given above, 

e. This study shows that gravity anomalies computed from levelling data 
related to either the Jervis Bay or Galveston Datums can bo assumed 
to refer to the geold to 0. 1 mGal, 

f. Tills technique also provides an exacting test of the value of GM used 
in Equation 40. 
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The Differential Effect of the Atmosphere in Geopotential Computations from Satellite 
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Table 2 

Statistics Froni Area Mean Values of (Fquation (16)) In the 
Australian Gravity Data Dank (AUSGAD 76), Based on the Freely 
Adjusted I.evel Network for Australia Referred to the Jervis Bay 
Datum Level Surface (Units mGal) 


d0 - 30* dx - 45* 


Square Size 

1* X 1* 

2* X 2* 

6* X 5* 

1 

No, of Blocks 

722 

181* 

30* (15) 

^ * 

1 Mean Value 

(Arj) 

3,65 

3.81 

3.41 (-0.18) 

mis 

16.6 

13.3 

8.3 (6.5) 

j Expected 1 f 1 

0. 02 

1 

0.02 

0. 04 (0. 05) 


♦Minimum Representation = 40 percent 


(Fitp.ires within brackets for 5* x 5* squares are based on a sample which 
Includes only squares w'here mean is computed from 25 1* x 1* values; 
l.e. , 100^ representation) 
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Table 3 


Statistics from Area Mean Values of Agj (Equation (16)) in the Gravity 
Data Set for Central North America, Based on Geopotential Estimates 
Related to Tlie Galveston Datum Level Surface (Units mGal) 

d0 » 20“ dx = 45“ 


Square Size 

1“ X 1“ 

2“ X 2" 

6“ X 5" 

No, of Blocks 

835 

218* 

34* 

Mean Value 
(^&d) 

-0. 76 

-0,67 

-0.25 

rms 

15,9 

10.4 

6.5 

Expected 1 J'sj / Ag j 1 

O.OS) 

0.03 

0.05 


*Minimum Representation = 40 percent 



















Table 4 

The Potential of the Geoid (Wq ) from GEOS-3 Altimeti'y 
GM » 3. 986 004 7 x lO^Ocm^ s-2 


Data Source 

Wallops 

Wallops 

Epoch 

Sep-Oct '75 
Mar-Apr '76 

Peb-Aug '76 

No, of Passes 

634 

882 

No, of 1® Sq, Sampled 

13,499 

12,349 

rms (Wjs - Wp ) A kGalm 

5.8 

5,1 

Wo(lcGalm) 

6,263,682,76 

6,263,682. 39 
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